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Abstract Very sensitive HI absorption spectra (At ~ 10"'' over 1 km/s) toward 
high latitude QSOs have revealed a population of tiny discrete features 
in the diffuse ISM with peak r of 0.1 - 2% and core line-widths corre- 
sponding to temperatures as low as 20 K. Imaging detections confirm 
I linear dimensions of a few 1000 AU. We suggest these structures may 

. be formed by the stellar winds of intermediate mass stars. A more 

' speculative origin might involve molecular "dark matter" . 

o 

^ ' 1. Introduction 

Over the past decades there have been several hnes of evidence sug- 
O I gesting a surprising degree of small-scale structure in the atomic inter- 

' stellar medium. One of the first of these was the observation of spa- 

. tially variable HI absorption seen toward compact radio sources (eg. 

Dieter et al. 1976, Davis et al. 1996, Faison et al. 1998, Faison & 
Goss 2001). In the most extreme cases, for example toward 3C138, 
there is good evidence for a variation in the HI opacity, r = 1.7(n • 
; s)/{TAV)[cm~^pc/K - km/s], of as much as At = 0.1 on 20 AU 

transverse scales (Faison et al. 1998). The "simplest" interpretation 
of these observations are very large variations in the volume density, 
AiiHi ~ 10^ cm~^, assuming that all of the other relevant variables 
(specifically pathlength, s, and temperature, T) are kept fixed. How- 
ever, as argued by Deshpande (2000), realistic ISM structure functions 
can lead to large variations of r with small angular offsets simply from 
the statistical fluctuations in the effective pathlength with position. An- 
other line of evidence for small-scale atomic structure came from searches 
for time variability in the HI absorption seen toward pulsars (eg. Frail 
et al. 1994, Johnston et al. 2003, Stanimirovic et al. 2003). However, 
the early claims for ubiquitous and significant time variations in r have 
not been confirmed in the recent careful studies. A third line of evidence 
has come from observations of Na I observation toward nearby pairs of 
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stars (Watson & Meyer 1996, Lauroesch et al. 1998). These observations 
have revealed strong evidence for discrete absorption features in the local 
ISM which have highly variable properties on lOO's of AU scales. Since 
discrete features are being probed by these observations, rather than 
simply the projection along a long line-of-sight, the structure function 
arguments of Deshpande do not seem to apply. 

2. Observations 

We have undertaken a series of extremely sensitive HI absorption ob- 
servations toward bright background QSOs near the North Galactic Pole 
(NGP) utilizing the Westerbork Synthesis Radio Telescope (WSRT). 
The initial motivation for these observations was detection of HI absorp- 
tion from a Warm Neutral Medium (WNM), even if it's temperature was 
as high as 10^ K. Such an experiment was prompted by the detection of 
HI absorption by Kanekar et al. (2003) at velocity widths extending at 
least as high as an equivalent temperature of 3500 K, which corresponded 
to the sensitivity limit of those data. The NGP region was chosen since 
it was expected that the lincs-of-sight through our Galaxy disk might be 
as short as possible and therefore relatively simple. Sensitive data were 
acquired toward 3C286 (14.7 Jy), 3C287 (7 Jy), 4C+32.44 (5 Jy) and 
B2 1325-1-32 (1.4 Jy) by observing in an in-band frequency-switching 
mode utilizing a 1 MHz throw every 5 minutes inside a 2.5 MHz to- 
tal bandwidth with 0.5 km/s channel width. The in-band frequency 
switching allowed exceptionally good band-pass calibration while pro- 
viding 100% of the observing time on-source. For the brightest sources 
we achieve an RMS Ar < 2 x 10~^ over only 1 km/s, making these 
the most sensitive HI (21cm) absorption measurements of which we are 
aware. 

3. Results 

Rather than detecting only the simple, broad and shallow absorption 
profile of the WNM toward these high latitude lines-of-sight, we were 
surprised to find instead that multiple narrow absorption features were 
detected at discrete linc-of-sight velocities with peak opacities in the 
range 0.1 to 2%. The detection of discrete absorption features was par- 
ticularly surprising since the peak emission brightness seen in these di- 
rections (with the 35 arcmin total power beam) was only between about 
2 and 5 K. The four observed lincs-of-sight are compared in Figure 1. 
Since many of the absorption lines are detected with high signal-to-noise 
it is possible to say with confidence that the profiles of individual features 
are non-Gaussian, and instead appear to be semi-Lorentzian with very 
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narrow line-cores that merge smoothly into broader wings. Of course 
it is always possible to represent each feature with a sum of coincident 
Gaussians, but such a procedure appears to be quite arbitrary and it is 
far from clear what physical interpretation might be attached to such 
an arbitrary decomposition. 
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Figure 1. HI absorption (left) and total power emission (right) spectra. 

Moderately near-by lines-of-sight show essentially uncorrelated ab- 
sorption features and only weakly correlated total power emission fea- 
tures. Only in the case of 4C-I-32.44 and B2 1325-1-32 which have an 
angular separation of only 15 arcmin is there some possibility of similar 
velocity components having been detected. 

4. Analysis 

Motivated by the non-Gaussian line profiles we have explored some 
simple spherically symmetric, iso-baric cloud models of the form: 

nH{r) = noexp[-{r/s)'^^] for T < 4000 K and (1) 

nnir) = noexp[-{r / s^'^] for T > 4000 K (2) 

where we further relate volume density to temperature using the ther- 
mal pressure which was assumed to be constant at Pjk^ = tihT = 
1500 cm~'^K. The temperature was allowed to vary between T^m = 
20 K and T^ax = 15000 K, yielding an assumed thermal velocity dis- 
persion of cr^ = 0.0086T. The predicted HI absorption and total power 
emission spectra were calculated for a "cloud" placed at the central ve- 
locity of each observed feature in an attempt to simultaneously reproduce 
the observed spectra shown in Fig. 1. The most important free parame- 
ters in this process were the cloud scale-length, s, the cloud distance, d 
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(which most strongly influences the predicted total power emission) and 
an impact parameter, b. This last parameter was used to allow for the 
likely circumstance that each spherical model cloud may not be pene- 
trated exactly on-axis by the background absorber. The two power-law 
indices of the scaled distance, ai and 0:2 determine the characteristic 
line-shape in the cold core and warm halo respectively. Although these 
were in principle also free parameters, it was found that only minor 
variations from "standard" values of about ai=l/4 and a2=l/8 were 
needed. 




Figure 2. Overlay of observed and modelled spectra for the 3C286 1-o-s (loft). 
Density and temperature profiles of all modelled "clouds" (right). 



An illustration of the simultaneous emission and absorption line fitting 
is shown in Fig. 2 for the 3C286 line-of-sight. No fitting has actually 
been carried out, but merely a x-by-eye to illustrate the possibilities of 
this approach. The entire set of density and temperature profiles for the 
required clouds along all four lines-of-sight is also shown in the Figure. 
This illustrates the basic similarities of the modelled features, although 
the apparent scale-lengths do vary by about an order of magnitude. The 
apparent distances used in reproducing the spectra was about 10 pc, 
and the typical half-density radius was 100 AU, although these were not 
very well-constrained given the very low angular resolution of our total 
power data (35 arcmin). These modelling results can be easily scaled 
to other assumed thermal pressures by a linear scaling of nn together 
with an inverse linear scaling of both s and d. For example, with a 
typical thermal pressure of only 150 cm~^K, comparable fits would be 
obtained at apparent distances of 100 pc and typical half-density radii of 
1000 AU. Variations from the idealized spherical cloud symmetry would 
have a similar impact on apparent distances and sizes. 
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We also have more concrete information regarding the nature of these 
features, since the same WSRT observations that provided the absorp- 
tion spectra also allow an imaging search for emission counterparts. 
Cubes of HI line emission were produced at a range of angular reso- 
lutions (15, 30, 60 and 120 arcsec). Obtaining sufficient brightness sen- 
sitivity to achieve detections in emission typically required smoothing 
to 60 arcsec. At those velocities where total power emission exceeding 
about 1 K is seen in Fig. 1 compact emission clumps of 2-3 K bright- 
ness are detected at apparently random locations in the field superposed 
on the poorly-sampled diffuse background emission. One such emis- 
sion clump, immediately adjacent to the 3C286 linc-of-sight is shown in 
Fig. 3. The intrinsic angular size of these clumps appears to be about 
30 arcsec, while their FWHM line- widths are 1-2 km/s, corresponding 
to the thermal lincwidths of 20-80 K HI. A representative detected col- 
umn density for the clumps in the 3C286 field is Nhj ~ 5 x lO^^cm"^. 
At an assumed distance of say, 100 pc, the clump size would correspond 
to 3000 AU and the central volume density would be about 100 cm~^. 

A more interesting emission structure is detected in the 4C+32.44 
field. The linc-of-sight toward this background source appears to in- 
tersect a 15 arcmin diameter shell of HI emission as shown in Fig.4. 
Although it may be a chance super-position, this apparent shell includes 
the GOIII star HD 116856 at (a, 5)2000 = (13:25:55.835,-^31:51:40.629). 
The measured parallax of this star places it at 105ibll pc, where the 
shell would have a 0.45 pc diameter. The stellar proper motion (Aa, Ad) 
= (-1-14.73,— 43.40) mas/yr is directed predominatly toward the South. 
Peak column densities in this structure reach Nhj > lO^^cm"^ with 
FWHM line- widths of 2-3 km/s. 
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Figure 4- HI emission shell toward 4C+32.44 in map (left) and spectrum (right). 



5. Discussion 

Although more work needs to be done to fully characterize the type 
and quantity of sub-structure in the "diffuse" ISM, it is already be- 
coming clear that even the most diffuse regions are populated by tiny 
distinct structures of very high density- and temperature-contrast. For 
example, the low scale-size end of the turbulent power spectrum is pre- 
dicted (Deshpande 2000) to have HI opacity fluctuations of only 10~^ on 
spatial scales of 1000 AU. In fact, we measure more than two orders of 
magnitude larger opacity fluctuations of > 10~^ on these scales. There 
appears to be substantial injection of fluctuation power on very small 
scales. The physical origin of these tiny structures is not yet clear, but 
it seems conceivable that the stellar winds of intermediate mass stars 
may play an important role in their formation, whenever such stars find 
themselves within a diffuse atomic structure. A more speculative origin 
might be some relation to molecular "dark matter" (eg. Pfenniger &; 
Combes 1994). New observations should test various scenarios. 
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